On the energy spectrum of strong magnetohydrodynamic turbulence

The energy spectrum of magnetohydrodynamic turbulence attracts interest due to its fundamental importance and its relevance for interpreting astrophysical data. Here we present measurements of the energy spectra from a series of high-resolution direct numerical simulations of MHD turbulence with a strong guide field and for increasing Reynolds number. The presented simulations, with numerical resolutions up to 2048^3 mesh points and statistics accumulated over 30 to 150 eddy turnover times, constitute, to the best of our knowledge, the largest statistical sample of steady state MHD turbulence to date. We study both the balanced case, where the energies associated with Alfv\'en modes propagating in opposite directions along the guide field, E^+ and $E^-, are equal, and the imbalanced case where the energies are different. In the balanced case, we find that the energy spectrum converges to a power law with exponent -3/2 as the Reynolds number is increased, consistent with phenomenological models that include scale-dependent dynamic alignment. For the imbalanced case, with E^+>E^-, the simulations show that E^- ~ k_{\perp}^{-3/2} for all Reynolds numbers considered, while E^+ has a slightly steeper spectrum at small Re. As the Reynolds number increases, E^+ flattens. Since both E^+ and E^- are pinned at the dissipation scale and anchored at the driving scales, we postulate that at sufficiently high Re the spectra will become parallel in the inertial range and scale as E^+ ~ E^- ~ k_{\perp}^{-3/2}. Questions regarding the universality of the spectrum and the value of the "Kolmogorov constant" are discussed.Comment: 13 pages, 10 figures, accepted for publication in Physical Review X (PRX

Scaling properties of small-scale fluctuations in magnetohydrodynamic turbulence

This is the final version of the article. Available from IOP via the DOI in this record.Magnetohydrodynamic (MHD) turbulence in the majority of natural systems, including the interstellar medium, the solar corona, and the solar wind, has Reynolds numbers far exceeding the Reynolds numbers achievable in numerical experiments. Much attention is therefore drawn to the universal scaling properties of small-scale fluctuations, which can be reliably measured in the simulations and then extrapolated to astrophysical scales. However, in contrast with hydrodynamic turbulence, where the universal structure of the inertial and dissipation intervals is described by the Kolmogorov self-similarity, the scaling for MHD turbulence cannot be established based solely on dimensional arguments due to the presence of an intrinsic velocity scale-the Alfvén velocity. In this Letter, we demonstrate that the Kolmogorov first self-similarity hypothesis cannot be formulated for MHD turbulence in the same way it is formulated for the hydrodynamic case. Besides profound consequences for the analytical consideration, this also imposes stringent conditions on numerical studies of MHD turbulence. In contrast with the hydrodynamic case, the discretization scale in numerical simulations of MHD turbulence should decrease faster than the dissipation scale, in order for the simulations to remain resolved as the Reynolds number increases

Dynamic Alignment and Exact Scaling Laws in MHD Turbulence

Magnetohydrodynamic (MHD) turbulence is pervasive in astrophysical systems. Recent high-resolution numerical simulations suggest that the energy spectrum of strong incompressible MHD turbulence is $E(k_{\perp})\propto k_{\perp}^{-3/2}$. So far, there has been no phenomenological theory that simultaneously explains this spectrum and satisfies the exact analytic relations for MHD turbulence due to Politano & Pouquet. Indeed, the Politano-Pouquet relations are often invoked to suggest that the spectrum of MHD turbulence instead has the Kolmogorov scaling -5/3. Using geometrical arguments and numerical tests, here we analyze this seeming contradiction and demonstrate that the -3/2 scaling and the Politano-Pouquet relations are reconciled by the phenomenon of scale-dependent dynamic alignment that was recently discovered in MHD turbulence.Comment: Published versio

Magnetic dynamo action in random flows with zero and finite correlation times

Hydromagnetic dynamo theory provides the prevailing theoretical description for the origin of magnetic fields in the universe. Here we consider the problem of kinematic, small-scale dynamo action driven by a random, incompressible, non-helical, homogeneous and isotropic flow. In the Kazantsev dynamo model the statistics of the driving flow are assumed to be instantaneously correlated in time. Here we compare the results of the model with the dynamo properties of a simulated flow that has equivalent spatial characteristics as the Kazantsev flow but different temporal statistics. In particular, the simulated flow is a solution of the forced Navier-Stokes equations and hence has a finite correlation time. We find that the Kazantsev model typically predicts a larger magnetic growth rate and a magnetic spectrum that peaks at smaller scales. However, we show that by filtering the diffusivity spectrum at small scales it is possible to bring the growth rates into agreement and simultaneously align the magnetic spectra

The electromotive force in multi-scale flows at high magnetic Reynolds number

Recent advances in dynamo theory have been made by examining the competition between small and large-scale dynamos at high magnetic Reynolds number Rm. Small-scale dynamos rely on the presence of chaotic stretching whilst the generation of large-scale fields occurs in flows lacking reflectional symmetry via a systematic electromotive force (emf). In this paper we discuss how the statistics of the emf (at high Rm) depend on the properties of the multi-scale velocity that is generating it. In particular, we determine that different scales of flow have different contributions to the statistics of the emf, with smaller-scales contributing to the mean without increasing the variance. Moreover we determine when scales in such a flow act independently in their contribution to the emf. We further examine the role of large-scale shear in modifying the emf. We conjecture that the distribution of the emf, and not simply the mean, largely determines the dominant scale of the magnetic field generated by the flow

Strong-field dynamo action in rapidly rotating convection with no inertia

The earth's magnetic field is generated by dynamo action driven by convection in the outer core. For numerical reasons, inertial and viscous forces play an important role in geodynamo models; however, the primary dynamical balance in the earth's core is believed to be between buoyancy, Coriolis, and magnetic forces. The hope has been that by setting the Ekman number to be as small as computationally feasible, an asymptotic regime would be reached in which the correct force balance is achieved. However, recent analyses of geodynamo models suggest that the desired balance has still not yet been attained. Here we adopt a complementary approach consisting of a model of rapidly rotating convection in which inertial forces are neglected from the outset. Within this framework we are able to construct a branch of solutions in which the dynamo generates a strong magnetic field that satisfies the expected force balance. The resulting strongly magnetized convection is dramatically different from the corresponding solutions in which the field is weak